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A STUDY OF THE DYNAMICS OF AIRPLANE BRAKING SYSTEMS 


AS AFFECTED BY TIRE ELASTICITY AND BRAKE RESPONSE 

By Sidney A, Batterson 
Langley Research Center 


SUMMARY 


An analog computer study of an automatic airplane hraking system vas made 
in which the effect of tire elasticity and the time required to apply and 
release brake torque were considered. The equations of motion were derived in 
general terms for a simplified single-wheel arrangement and then solved on the 
con 5 )uter for the particular case of a wheel equipped with a 52 x 8.8 type VTI 
tire. It was found that the braking-system efficiency depended on a complex 
interaction between the values of tire frequency^ brake response time, and skid- 
control-sensor frequency, and depending on the combinations of these values, 
the braking -system efficiency could vary from 0 to 90 percent. The solutions 
also Indicated that a braking system which operated satisfactorily in the pres- 
ence of relatively hi^ or medium coefficients of friction could operate unsat- 
isfactorily at low coefficients of friction. It was also found that for low 
friction operation, large decreases in both stopping distance and tire skidding 
could be realized by proper selection of the wheel deceleration value used to 
generate the brake -release signal. 


INTRODUCTION 


The braking systems of most high-performance aircraft in use today utilize 
automatic anti-skid, or more appropriately, skid-control devices which regulate 
the torque applied by the brakes during the landing run. On dry and damp run- 
ways where the available coefficient of friction is large, the combined anti- 
skid control and brake system is efficient in minimizing airplane stopping dis- 
tance and also is effective in preventing excessive tire skidding. On water- 
and slush-covered runways, however, where the available coefficient of friction 
can be small, a number of hazardous landing Incidents have occurred which were 
attributed to inadequate braking. (See, for example, ref. 1.) 

In an effort to obtain further Insight into the phenomena which affect the 
basic operation of skid-control systems and also to determine the effect of low 
tire-ground friction coefficients on these systems during braking, an analog 
computer study was carried out for a simplified mathematical model of a braking 
system consisting of a skid control, brake, wheel, and tire. In this study, an 
effort was made to account for both the elastic behavior of the tire and the 
brake time response during cyclic braking. 



With regard to the mechanical simulation programed for the computer,, it 
should he pointed out that very little is known about the general dynamic behav- 
ior of a rolling wheel and elastic tire subjected to brake torque. This is par- 
ticularly true of the mechanical properties of the tire and the forces that 
occur at the tire-ground interface. Therefore, it was not possible in this 
study to represent analytically some of the effects which are known to occur 
during braking; for example, the change in tire-rolling radius caused by 
stretching of the rubber, the variable location of the vertical center of gro\ind 
pressure, and the development of ground force in the tire footprint area. 
Nevertheless, it is shown that the simulation chosen for this study resembles 
the observed dynamic behavior of actual wheels and tires subjected to braking 
forces . 

The mathematical model of the braking system consisted of a rigid wheel 
and rigid tire connected by a spring and viscous damper to simulate the elastic 
eind hysteresis properties of a pneumatic tire. The brake torque was programed 
to Increase and decrease with time at various specified rates, and the signal to 
apply and release the brake torque was generated by a control that sensed wheel 
angular acceleration. 

The basic dynamic equations for the braking system were derived in general 
terms and then solved on the analog computer for the specific case of a wheel 
and a 52 X 8.8 type VII tire having an initial horizontal velocity of 200 feet 
per second (6O.96 m/s). The tire radius used in the computations was I.I6 feet 
(0.5556 m), and a constant vertical ground reaction of 22,000 pounds (97,856 N) 
was applied to the tire throughout the entire landing run. Solutions were made 
for a range of tire-runway friction coefficients, brake time responses, control 
frequencies, and acceleration values at which the control initiated the applica- 
tion and release of brake torque. The results indicated the effect of these 
parameters on stopping distance and tire skidding. 


SYMBOLS 


Measurements for this investigation were taken in the U.S. Customary System 
of Units. Equivalent values are indicated herein in the International System 
(si) in the interest of promoting use of this system in future NASA reports. 
Details concerning the use of SI, together with physical constants and conver- 
sion factors, are given in reference 2. 
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/ w \ 
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drag load tetween tire and ground^ Ibf (n) 
acceleration due to gravity^ ft/sec2 (m/sec^) 
moment of inertia of wheel ^ ft-lbf-sec^ (kg-m^) 

moment of inertia of tire, ft-lbf-sec^ (kg-m^) 

moment of inertia of control sensor, ft-lbf-sec^ (kg-m^) 

spring constant of tire, Ibf /ft (N/m) 

spring constant of control sensor, Ibf/ft (N/m) 

control sensor mass radius, ft (m) 
generalized force where i = 0, cp, and p 

generalized coordinate where i = x, 0, cp, and P 
wheel radius, ft (m) 
tire radius, ft (m) 

effective tire-rolling radius during braking, ft (m) 
effective tire-rolling radius during free rolling, ft (m) 
slip ratio of wheel 
slip ratio of tire 

total kinetic energy of braking system, ft-lbf (m-N) 
braking torque, ft-lbf (m-N) 
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braking torque during wheel locked condition, ft-lbf (m-N) 
maxlmiim braking torque, ft-lbf (m-N) 
time, sec 

time for brake torque to Increase from zero to Its maximum, sec 
time for brake torque to decrease from its maximum value to zero, sec 
potential energy, ft-lbf (m-N) 

potential energy of tire due to elastic deformation of tire , ft-lbf 
(m-N) 

potential energy of sensor mass due to elastic deformation of sensor 
spring, ft-lbf (m-N) 

total wei^t, Ibf (N) 

work done by tire in skidding, ft-lbf (m-N) 

linear displacement of wheel axle, ft (m) 
initial forward velocity, ft/sec (m/sec) 
tire circular spring deflection, ft (m) 
control sensor circular spring deflection, ft (m) 


angular displacement of control sensor mass, rad 

control sensor acceleration value that generates brake application 
signal, rad/sec2 

control sensor acceleration value that generates brake release sig- 
nal, rad/sec2 


braking-system efficiency. 


Minimum stopping distance ^ percent 

Actioal stopping distance 


angular displacement of wheel, rad 
coefficient of friction 

maximum obtainable coefficient of friction between tire and ground 
coefficient of friction at zero tire angular velocity 


k 


^0 

Icp 


torsional damping ratio of wheel due to tire 
torsional damping ratio of tire 
torsional damping ratio of control sensor 


skid index, 


Total work done by tire in skidding ^ percent 

Initial airplane kinetic energy 
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angular displacement of tire^ rad 
natural frequency of wheel on tire ^ rad/sec 

natural frequency of tire on wheel, rad/ sec 
natural frequency of sensor mass, rad/sec 


Dots over symbols indicate differentiation with respect to time. 


MATHEMATICAL MODEL AND METHOD OF ANALYSIS 


The physical system programed for the analog computer is shown schematically 
in flgiJire 1. It consists of a rolling wheel and tire with the total mass con- 
centrated at the axle and is traveling in the x-direction. The wheel is rigid 
and for the purpose of this analysis the tire is also assumed to he rigid and is 
connected to the wheel by a linear spring and dashpot in order to represent the 
elastic and damping character- 


istics of the tire. The wheel 
is subjected to a braking 
torque Tg which causes rela- 
tive motion between the wheel 
and tire. The resulting 
spring, damping, and tire 
inertia forces produced by 
this motion causes a drag 
force D to act between. the 
tire and the ground. This 
drag force is expressed as 
D = |iW where \i is the coef- 
ficient of friction between 
the tire and the ground, and 
W is the total weight. The 
angular displacement of the 
wheel is designated 0 and 
the angular displacement of 



the tire is designated q). 
The skid-control-sensor unit 


Figure 1.- Mathematical model of braking system. 
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is siinuia.tc(i s- "balsinccd cylindrical mass free io pivot- about thic axlo center 
line and is restrained by a spring and dashpot attached to the wheel as shown 
schematically in figure 1, The angular displacement of the sensor mass is 
designated p. The computer was programed to generate brake application and 
brake release signals at certain preselected values of P (the angular accel- 
eration of the sensor mass). 

Althou^ current skid-control systems are much more sophisticated than this 
acceleration type^ it was selected on the basis of simplicity for this first 
analysis. The operation of this type of control unit is identical to the opera- 
tion of many standard angular accelerometer Instruments and some of the first 
on-off skid-control units operated in this manner. It is assumed that through- 
out the braking run, the vertical ground reaction is directly beneath the wheel 
center and remains equal to the total weight. 


Equations of Motion 

The variables x, 0, cp, and p are considered as generalized coordi- 
nates, and it is shown in appendix A that the equations of motion of the system 
are 


X = -pg 


T - 0) + - 9) -- 

3 = -203n,p^(3 - 0) - ^,p(3 - 0) 



(1) 

= Zs 

(2) 

% 

M-^^T 

(5) 

i-r 

0) 

w 


In order to obtain a solution for these equations, it is necessary to 
define a set of initial conditions and also to define a variation for q a^nd 
Tg in terms of known quantities c The initial conditions are obtained by 
selecting an initial forward speed Xq which then defines the initial values of 
0, cp, and p through the kinematic relation for free rolling. The initial 
displacements x, 6, cp, and p are taken equal to zero. 


Coefficient of Friction 

The variation of the coefficient of friction is defined in a somewhat dif- 
ferent manner than is followed in current practice. As pointed out previously, 
the wheel and the tire can have different motions because the tire is considered 
to deform elastically under the application of brake torque (simulated in this 
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case by the spring connecting the tire and wheel), whereas the wheel is assumed 
to be rigid. It is clear that since the coefficient of friction is a function 
of the relative motion between two rubbing surfaces, the coefficient of friction 
between the tire and runway must be a function of the tire motion cp with 
respect to the runway rather than of the wheel motion 9. It is customary in 
current practice to define the variation of the coefficient of friction pi as 
a function of the slip ratio, where the slip ratio is defined as 


% = 



( 5 ) 


where Sy is the slip ratio defined as a fxinction of wheel angular velocity, 9 
is the wheel angular velocity, x is the axle linear velocity, and Tq is the 
tire-rolling radius . 


The wheel angular velocity is used in equation (5) since in experimental 
testing this quantity can be easily measured during the application of brake 
torque, whereas the velocity of the tire in the footprint would be very diffi- 
cult to obtain. The equations of motion (eqs. (l) to ( 4 )), however, provide a 
solution for the tire velocity cp; therefore, in this analysis the slip ratio 
between the tire and the ground is used and is defined in a somewhat similar 
manner as 


Siji 


1 


y 

x/nj. 


( 6 ) 


where is the slip ratio defined as a function of tire angular velocity^ (p 

is the tire angular velocity^ and rrp is the tire radius. (See fig. 1.) 

In order to define the relation between |i and it was assumed that 

as torque was applied to the freely rolling wheel , no slip occurred between the 
tire and runway surface until the maximum apparent coefficient of friction was 
developed. The maximum apparent coefficient is defined as the maximum coeffi- 
cient of friction that could be developed at the tire -ground interface at the 
Instant under consideration. This maximum coefficient is indicated symbolically 
as f-imax* assumption of no slip (that is, Srp = O) from free rolling until 

the instant when is reached presupposes that the coefficient developed 

within this range depends primarily on tire stiffness, and this supposition is 
in general agreement with the opinions expressed by a number of investigators. 
(See refs. 5, 4, and 5.) Some experimental results, which also tend to bear 
out this simulation, will be discussed subsequently. 

If the brake torque is increased following the attainment of M^nax^ 
tire starts slipping with respect to the ground and the slip ratio Sqi takes 
on increasing values. A value of Sip of 1 indicates that the tire is in a full 
skid and has no rotational motion j that is, 9 = 0. The coefficient of friction 
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associated with a full skid is indicated symbolically as general, 

for the conditions of interest in the stopping of airplanes, \^jnax con- 

sidered greater than 4skid* 


The manner in which the coefficient of friction varies with slip ratio for 
coefficients lying in the region between and depends on a large 

number of variables which include, for exairple, surface roughness, tire tread, 
tire and ground temperature^ runway contamination, tire pressure, and so forth- 
Because little is known of the manner in which the coefficient varies with 
many of these parameters, and because some vaiy in a random manner throughout 
a braking stop, an empirical relation was formulated for use in the analog- 
computer study. To obtain this relation, use was made of experimental braking 
data obtained during the tests reported in reference 6. The data were normal- 
ized by using the ratio of instantaneous friction coefficient to maximum coef- 
ficient h/hmax" experimental data points which were obtained during the 

first braking cycle of some of the dry surface runs are plotted in figure 2. 

As would be expected, these data scatter over a wide band because, as pointed 



Figure 2.- Variation of friction coeffi- 
cient vith wheel slip ratio obtained 
from experimental braJ^ing test on a 
dry concrete surface with a 52 X 8.8 
type VII tire, (Data obtained from 
ref , 6 . ) 

occurs (sometimes referred to as the 
tire) that is, for values of S«p > 0 


out previously, the value of the coef- 
ficient depends on a number of vari- 
ables, some of which vary in an unpre- 
dictable manner. For this reason 
there is no unique curve from a prac- 
tical standpoint that describes the 
variation of p with Sip for values 
of Sp > 0. The variation for Sp > 0 
chosen for the analog computations was, 
therefore, obtained by choosing an 
expression which gave an average 
fairing of the experimental data. 

The variation of p/pmax 
Sp used in the calculations is shown 
in figure 5(a). Since Sp is the 
slip ratio referred to the tire, and 
since it was assumed that no slip 
occurred in the tire footprint until 
p reached the value of Pmax^ 
be seen that the slip Sp is equal to 
zero in this region (from free rolling 
where p = 0 until p = p^^^^). The 

empirical expression used to define the 
variation to the right of where Pjj^^x 

ck side of the p-slip curve for the 


M-max 



( 7 ) 
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(a) Analog simulation of variation of (t) Variation of friction coefficient 

friction coefficient with tire slip ratio. with wheel slip ratio. 

Figure 5.- Experimental hraking data and analog outputs of the variation of fric- 
tion coefficient with tire and wheel slip ratios for an Initial braking cycle. 

For values of Sip > 1 the value of the coefficient of friction was taken equal 
to that for S<p = 1. 


Figure 3(h) is an actual curve obtained from the analog computer of a 
first-cycle braking variation using the simulation shown in figure 1. In fig- 
ure 3(b) M-/|dmax plotted against Sy, the slip ratio referred to the wheels 
which is the manner in which data on the variation of with slip ratio are 
commonly presented. It can be seen that this initial braking cycle has^ in 
general^ the form of curves usually obtained from experimental braking tests; 
that is, Mmax occurs at a slip ratio greater than zero and the values of Sy 

to the left of the peak (sometimes referred to as the front side of the |^-slip 
curve for the wheel) depend primarily on the tire stiffness. 


To justify further the simulation of fi with S«p used for this analog 
study in the region between |i = 0 and u - 4max^ some data from the experi- 
mental braking tests reported in reference 6 are presented in figure k. The 
curve associated with the square symbols was obtained during a free-rolling 
test at a forward speed of about I 50 ft/ sec (i+ 5.72 m/sec); that is, no braking 
torque was applied to the wheel. The other curve was obtained for a test made 
at the same forward speed by using the same 32 X 8.8 type VII tire; however, 
during this test, brake torque was being applied to the wheel and was 
increasing approximately linearly with time throu^out the distance indicated 
by the data. At the start of the measurements, the brake-torque value was 
about one-quarter of that required for locking the wheel, whereas at the time 
that the last data point was obtained, the brake torque was approaching the 
value necessary to develop M-max’ instrumentation used for these tests 

provided for extreme accuracy in measuring and correlating the parameters of 
distance traveled x and wheel angular displacement 0. The static vertical 
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Distance along runwayj x, ft 

_1 ] J 

30 60 90 

Distance along runway, x, tn 


Flgiire Experimental variation of wheel rotation with dis- 
tance traveled obtained diiring a braking test and a free- 
rolling test made with a 32 X 8.8 type VII tire. Static 
load on tire^ 10,000 lb (^4,800 N); horizontal velocity, 

150 ft/sec ( 45*72 m/sec). (Data obtained from tests of 
ref. 6.) 

load for both tests was 10^000 pounds ( 44^480 N) and the tire pressure was 
260 pounds per square inch (179*26 N/cm^), it can be seen that the value of 
despite the increasing brake torque^ is apparently constant since the 
curves are straight lines and have a common origin. The value of Sy for the 
braking run can be computed from the equation given in reference 7 3-s 


Sw = 


- 


"e ,0 


^e 


where r^ is the effective tire-rolling radius during braking^ ahd r^ q 
the effective tire -rolling radius during free rolling. 


is 



The data presented in figure k yield values of of 1.^1 ft (0.5995 '^) 

and 0 of 1.25 ft (O. 38 IO m)j therefore, Sy = 0.046. 

In figure 5(t>) the computed value of Sy was also approximately constant 
for the upper 75 percent of the front side of the curve of variation of p 
with Sy (hereafter designated p-slip curve for the wheel) . Also the computed 
value for Sy in figure 3(^) is seen to average slightly over 0.05- Tho dotted 
line in figure 5(h) is the experimental hraking data of figure 4. This excel- 
lent agreement between the experimental braking results, and the results 
obtained from the analog computer would indicate that the variation of p with 
Srji chosen herein (fig. 3(a)) to simulate the action of an elastic tire during 
the application of brake torque gives a very good description of the actual 
variation. 

In order to use the variation of p with S<p shown in figure 3(a), it 
was necessary to program the computer to calculate p by two different methods. 
The particular method used at the instant the calculation was being made 
depended on whether Sip was equal to zero or greater than zero. During the 
time Sp = 0, the variable cp and its time derivatives were eliminated from 
equations (2) and ( 3 ) by using the expression 



which expresses the kinematic relation between the variables 9 and x when 
the tire is rolling without slipping, that is, Sp = 0. During this time, the 
computer operated in a mode which calculated the variables p, x, 0 , and p, 
and p for this case was a fiinction of tire elasticity, damping, and inertia. 
While in this mode of operation, the computer was programed to compare the 
value calculated for p with the value of Pmax given for the particular run. 
When the calculated value of p became greater than p^^ ^ , the computer 

switched over to a second mode of operation. In this second mode, the computer 
calculated the value of p as given by equation ( 7 ) and also calculated the 
values of x, 0, cp, and 3 by using equations (l) to (4). While in this 
mode of operation, the computer was programed to examine the value of Sp and 
continued to operate in this mode as long as Sp was greater than zero. When 
Sp became zero, the computer switched back to the first mode of operation. 

The computer then operated throughout the entire run and switched modes 
according to this described logic system. It should be pointed out that the 
value of p calculated in the first mode is not in a true sense a coefficient 
of friction since it does vary and no sliding occurs. It can, however, be 
interpreted as an apparent coefficient of friction. 


Brake Torque 

In this study, the brake was considered to be torque limited and the lim- 
iting value is indicated as further assumed that when brake 
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torque was applied or released, it varied linearly with time. Although the 
linear variation is never obtained in practice, the actual variation between 
brakes and brake installations can vary so widely that this linear variation 
was chosen in the interests of simplicity for this first analysis. The brake- 
torque — time variation used is illustrated schematically by the solid line in 
figure 5> The time for torque to increase from 0 to its limited value Tg 

is called tp and the time to decrease from Tg^j^g^ to 0 is called t2- Val- 
ues of tp and t2 could be selected at will in order to simulate brakes 
having various response times. 



Figure 5.- Variation of brake torque with time used In analog com- 

puter braking study. 


The computer was programed to permit the brakes to be released prior to 
reaching T3 they could be reapplied during the release interval before 

reaching zero torque. For example, refer to figure 5 3-iid suppose that the brake 
torque is increasing and at time tg, a release signal is generated by the 
accelerometer j at that instant the brake torque will start decreasing at a rate 
equal to that indicated by line 3—4-. Of course, if the release signal occurred 
some time after point 2 on line 1-2 was reached, the brake torque would con- 
tinue at this constant limiting value from the point 2 until the release signal 
was generated. If a brake application signal is generated before the torque 
reaches zero, for example, at time t^^, the torque will start increasing at a 
rate equal to that indicated by line 1-2. This cycling process is repeated by 
the computer throu^out the entire landing runout, the brake torque being sub- 
ject to the command signals generated by the control. 


Wheel Locking 

In many instances during the computer runs, the dynamic conditions of the 
problem were such as to allow the wheel to lock during certain portions of the 
runout j that is, the response of the brakes and/or the control was not fast 
enough in reducing brake torque to prevent the wheel from being braked to a 
complete stop. In order to eliminate prolonged periods of wheel locking, the 
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computer was programed to generate a "brake release signal when the wheel angu- 
lar velocity was less than 10 radians per second. Even so, some wheel locking 
did occur, and during this time the brake torque applied to the wheel was a 
function of the tire inertia and the applied ground drag torque. The locked- 
wheel brake torque is shown in appendix B to be 


'^B,L “ “ l^WrQi 


(8) 


When wheel locking occurred^ the computer was programed to compute hoth Tg 
and compare the two values. As long as T;g was greater than 

T;B^L> wheel remained locked and the brake torque used in the equations was 
Tb^L* When^ however^ Tb became less than Tb^L^ value of Tb vas used 

in the equations and the wheel began turning once again. During wheel-locked 
conditions^ ^-PP^strs as an oscillation in brake torque at the tire fre- 

quency because of the oscillatory nature of 9 immediately following wheel 
locking. (See fig. 6 .) 


Tire Skidding 

An important factor used for judging the performance of a braking system 
is the amount of tire skidding that occurs during the braking run. Tire wear, 
particularly at the higher coefficients of friction, depends largely on the 
amount of work done in skidding. Also, the amount of cornering force that can 
be developed is affected by tire skidding. The total work done in skidding is 
shown in appendix C to be: 


^skid ^ 



jJtSrjiX dt 


( 9 ) 


where 

Wgkid ^ork done in skidding 

W vertical wheel load 

[i instantaneous coefficient of friction 

Siji Instantaneous slip ratio referred to tire 

X instantaneous forward speed 

t time 


The analog computer was programed to solve equation (9) and produced a time 
history of the work done in skidding throughout the entire braking stop. 
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1 sec 


1 sec 



("b) Enlarged section from figure 6(a). 

Figure 6.- Concluded. 

COMPUTATIONAL PROCEDURE 
Physical Constants 

As indicated previously^ the analog computer vas programed to solve the 
dynamic equations of the braking system for a specific wheel, tire, and brake 
combination. The initial condition of forward speed of 200 feet per second 
(60.96 m/sec) was the same for all mins. Physical constants were chosen to 
describe a J2 X 8.8 type VII 22 ply rating tire and wheel combination. These 
constants are given in table I. This particular wheel and tire size was used 
because a large amount of braking data has been obtained for this combination 
at the Langley landing-loads track (ref. 6). It should be pointed out that the 
value of 1.16 feet (O.3556 m) chosen for r,p was an average of the unbraked 
tire-rolling radius and the deflected radius. For this tire, the unbraked 
rolling radius for a vertical load of 22,(X)0 pounds (97^856 N) is 1.22 feet 
(0.5719 whereas the deflected radius at this load is 1.10 feet (0.5553 • 

The value of the torsional spring constant of the tire krp was obtained from 
the empirical formula presented in reference 7* 
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TABLE I.- BRAKING-SYSTEM CONSTANTS 


Tire radius ; rrp, ft (m) . 

Tire moment of inertia , ft-lbf-sec^ (kg-m^) , . , 
Tire spring constant, kji, Ibf/ft (N/m) ....... 

Tire damping ratio, 

Wheel radius, ry, ft (m) 

Wheel moment of inertia, I^, ft-lhf-sec^ (kg-m^) . . 

Maximum brake torque, ft-lhf (m-N) 

Vertical load, W, Ihf (N) 

Initial horizontal velocity, Xq, ft/sec (m/sec) . . . 
Sensor control damping ratio, 


1.16 

(0.3536) 

1.875 

(2.542) 

120,000 

(175,080,000) 

0.1 


0.715 

(0.2175) 

0.855 

(1.1592) 

15,400 

(18,167.72) 

20,000 

(97.856) 

200 

(60.96) 

0.7 



Test Variables 


The test program consisted of obtaining computer solutions for some 
selected values of the brake time response tq and t2, the control accelera- 

9m ■ • 

tion at which brakes were applied and released (3q and angular 

accelerometer control frequency Certain combinations of these variables 

were used with the maximum coefficient of friction u 


max 


between tire and run- 


way varying over a range from 0.1 to 0.8. The control accelerometer damping 
coefficient |p was 0.7 for all runs.. The range of the brake and control 


parameters investigated was restricted because the number of these variables 
and the possible combinations of these variables, even though confined to a 
practical range, would involve a formidable con 5 )uter program if expanded to any 
extent. Because of the preliminary nature of this study, it was felt advisable 
to limit the scope of the investigation to a minimum number of solutions. 


Accuracy Checks 

In order to eliminate errors arising from unintentional circuitry changes 
and computer malfunctions, a series of frequency-response checks were carried 
out on individual components of the program each day, and also as a final check 
a standard test condition was run through the computer prior to each day*s oper- 
ation. The standard test condition incorporated the initial conditions and con- 
stants presented in table I. The other parameters associated with this run 
were 

^^max ^ ^ ^ 
tq = 0.1 sec 

t 2 = 0.1 sec 
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2jt 


= 0.5 cps 


32 = -50 rad/sec^ 

3']_ - 0 rad/sec^ 

The time-history traces of the computer outputs obtained for the standard run 
are shown in figure 6(a), Use of the standard run in checking consisted of a 
visual comparison of each channel with the corresponding channel of a master 
standard run. In addition^ the values obtained for the total runout x and 
the total work done in skidding were found to check the master values 

within about 2 percent when everything was operating properly. 

The saw-toothed appearance of x, the forward displacement time history^ 
resulted from the method programed to obtain a maximum output sensitivity from 
this channel. For this channel, the maximum galvanometer displacement was set 
to represent 100 feet (30,48 m) of displacement, and as each 100-foot (30.48-m) 
segment was traversed, the output polarity was reversed so that the saw-tooth 
record was obtained. Total displacement is obtained by multiplying the number 
of lines connecting the upper and lower peaks (in this case 20) by 100 and then 
adding the proportional deflection of the last partial-line segment. For the 
record shown in figure 6(a), the total runout was 2,075 feet (632.46 m) - 

For these runs the computer was slowed down so that its actual speed of 
operation was l/20th of real-time operation speed. This time scaling was nec- 
essary because the frequency response demanded of the computer in solving this 
problem was greater than its capability in real time. 

During some of the runs, selected outputs were recorded by a mechanical 
X,Y plotter. In order to Insure that the inertias of the X,Y-plotter components 
were not attenuating the plotted results, the plots were compared with observa- 
tions made of the identical run using a cathode ray oscilloscope. 


Calculation of Braking-System Efficiency and Skid Index 

Some results of the computer program are discussed in terms of braking- 
system efficiency and also in terms of the braking- system skid index. Effi- 
ciency r\ is defined herein as the ratio of the minimum stopping distance to 
the actual stopping distance multiplied by 100 to express it in percent. The 
minimum stopping distance depends on the initial velocity and the value of 
^max* actual stopping distance is the final value of x obtained from the 

computer. Expressed symbolically, the efficiency is 


T] = 




2 *^maxS 


X 


X 100 
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Skid index cr is defined as the ratio of the total work done in skidding 
to the initial kinetic energy of the airplane multiplied by 100 to express it 
in percent. Expressed symbolically, the skid index is 


2g 

The work done in skidding can be obtained from the computer results by using 
equation (9)^ and the initial kinetic energy depends on the initial airplane 
velocity and weight. This parameter is a measure of the amount of total 
stopping energy contributed by skidding of the tire. It then follows that for 
optimum operation of braking systems^ the efficiency should be as large as 
possible and the skid index as small as possible. 


RESULTS MD DISCUSSION 


The results and discussion are presented in two parts: The first part 

deals with the effects of tire elasticity^ brake time response ^ and control 
frequency on the braking system when the available tire-runway friction coef- 
ficient is large enough to be comparable with the coefficients encountered 
during landing stops on most dry or damp runways. The magnitude of the coeffi- 
cient used^ however^ is low enough to permit the maximum brake torque to lock 
the wheel; that is^ ‘^B^max ^ l^max^^T* second part discusses the effect of 

lower friction coefficients on the braking systems. Low coefficients may be 
experienced when runway surfaces are contaminated with foreign substances such 
as water ^ slushy snow^ and ice. 


Large Tire-Runway Friction Coefficients 

Effect of tire elasticity on wheel motion .- As was mentioned earlier when 
tire torsional elasticity is considered^ the wheel motion may differ from the 
motion of the tire in the ground-contact area. This condition is particularly 
true when brake torque is being cycled at a rapid rate. The significance of 
this differential motion between the wheel and tire lies in the fact that the 
skid-control sensor is mounted on the wheel whereas the vehicle retarding force 
is developed in the tire footprint and as pointed out previously, this force 
depends on the tire motion. Therefore, the control- and brake -response fre- 
quencies must be chosen so that the applied brake torque controls the tire 
motion in the desired manner even though the brake-control signals depend on 
the wheel motion. 

An example of the character of the individual motions of the wheel and the 
tire during cyclic braking can be seen in the standard run record presented in 
figure 6. During this ran, brake torque was cycled at a fairly rapid rate by 
the control and complete wheel locking begins after a runout of about 1,000 feet 
(50^,80 m) . In this region the wheel angular velocity 9 reaches zero and the 
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wheel remains locked for a brief period of time, (approximately 0.2 second); 
then as the brake torque drops off, it rapidly returns to the free-rolling con- 
dition. This condition of wheel locking followed by free rolling then continues 
in a cyclic manner throughout the remainder of the run. During the first 
1,000 feet {y^k.QO m) of the runout, however, the wheel angular velocity was 
large enou gh to make the response of the accelerometer control and brake ade- 
quate to prevent complete locking. Although complete locking did not occur in 
the first 1,000 feet (504.80 m) , it can be seen that both the wheel and tire 
slip-ratio traces indicated a significant amount of wheel and tire slip. 

The character of the relative motion between the wheel and tire can be 
observed by comparing the wheel and tire angular velocity traces 6 and cp 
in figure 6 during the period of wheel locking and spin-up. It can be seen 
from the 9 trace that the tire angular velocity becomes oscillatory and even 
exhibits negative values (rearward rotation) during some instants while the 
wheel is fully locked, that is, 9 = 0; or in other words, after the wheel came 
to a complete stop, the tire made several torsional oscillations about the axle 
center line. The oscillation in 9 while 0 has a zero value is a result of 
tire elasticity and tire inertia and can be seen on the tire slip-ratio trace 
Srp as an oscillatory variation which takes on values greater than 1 when the 
tire angular velocity is negative. The foregoing results indicate differences 
between wheel and tire motion when the wheel is locked and wheel and tire slip 
ratios have values in the neighborhood of 1. 

Differences between wheel and tire motion due to the elastic behavior of 
the tire can also be seen at the other end of the slip-ratio scale, that is, 
when Sy and 89 are in the neighborhood of zero. For this condition, how- 
ever, the wheel slip -ratio trace indicates that the slip ratio of the 

wheel takes on negative values, whereas the tire slip ratio is never less then 
zero. This result indicates that during wheel spln-up, the spring force devel- 
oped by the tire drove the wheel to angular velocities greater than ^ those cor- 
responding to that required for free rolling at the forward speed x. This 
phenomena of wheel angular velocity having larger values just after spin-up 
than it had prior to wheel locking is often observed during both landing-loads 
track and dynamometer -braking tests. It should be stated, however, that during 
actual braking tests, tire stretch occurring dxiring the application of drag 
load also adds to this effect. 

The overall effect of differences in wheel and tire motion on the varia- 
tion of n with slip ratio throu^out the entire slip-ratio range can be seen 
in figure 7. These curves (fig. 7) were obtained during the braking stop shown 
in figure 6 by connecting the computer outputs of friction coefficient and slip 
ratio to a mechanical X,Y plotter where [i was plotted against slip ratio. 

The plot on the right shows the programed variation of p with 89 and 
remained unchanged throughout the entire landing run; the plot on the left, 
however, which indicates the variation of p with 8^, exhibits significant 
variations from cycle to cycle. It is of particular significance to note that 
the slope of the curve of p plotted against 8w can have either positive or 
negative values when the value of 89 = 0. This result is in contrast to cur- 
rent notions regarding this variation since experimental data for p plotted 
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Slip ratio referred to wheel, 


Slip ratio referred to tire, Sj. 


Figure 7.- Analog computer outputs obtained over a large number of braking cycles showing variation 
of friction coefficient with wheel and tire slip ratios. Data obtained from test shown in fig- 
ure 6. 


against slip is normally presented with this slope drawn positive. This result 
would indicate that for the conditions of this run, a control sensor located 
on the wheel could not he used to determine whether the tire is operating in 
the range between free rolling and or between locked- 

wheel condition, just on the basis of the slope of the }i-slip curve for the 
wheel. 


Effect of tire elasticity on coefficient of friction .- Another aspect of 
tire elasticity is indicated in the coefficient of friction trace appearing in 
figure 6. During the time that the wheel is locked, the coefficient of fric- 
tion is substantially constant and equal to ^skid for a small spike 

caused by the tire oscillation. When the brake torque drops to a value which 
is low enough to permit the wheel to start turning, the value of [x rises very 
rapidly. This rapid rise is caused by the increase in the spin-up drag load as 
the tire goes from a full-slip to a zero-slip condition during the time that 
the brake torque is decreasing; that is, at the instant the wheel starts 
turning, the spin-up drag load is equal to 3<nd as the brake torque 

continues to decrease, the drag load increases until its value is M^max^ 
the instant the tire slip ratio Srp reaches zero. It is of interest to note 

that for the conditions of this run, the time period required for the wheel to 
spin up from its fully locked condition to free rolling was somewhat less than 
the time period required to go from the initial incipient skid condition to a 
fully locked wheel condition. This fact can be observed by coir^aring the slopes 
of both the \i and 0 time -hi story curves during spin-up and spin-down. 
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Immediately upon reaching hjnax after spin-up, it can he seen that the 

value of n drops very rapidly toward zero. The decrease is so rapid that the 
record trace becomes much less pronounced because of the hi^er galvanometer 
writing speed. The time rate at which n decreases in this run depends almost 
entirely on the natural period of the tire. It is evident from figure 6(b) 
that the coefficient of friction did not return to a zero value prior to the 
start of the next cycle. The reason for this effect was that even though this 
was a short period of time, it was long enough for this control and brake to 
reapply a sufficient amount of brake torque to prevent the friction coefficient 
from reaching a value of zero. 

The fact that the value of the coefficient of friction is affected by the 
response period of the tire indicates that the stopping distances would also be 
affected by tire elasticity. From a practical standpoint, this result would 
indicate that the same brake and control system could exhibit different stopping 
efficiencies when used with wheel and tire combinations having different tor- 
sional frequencies. 

Control-sensor and brake -time response . - Figure 8 shows the effect of var- 
iation of the accelerometer-control frequency on braking-system efficiency for 



Natural frequency of control sensor, cps 


Figure 8.- Variation of tra king- system efficiency with sensor- 

control frequency for three brakes having different torque decay 
rates. 3^^ = 0; ^2 ~ "30 radians/second^; =^0.4. 

three brakes By and C which differ from each other only in their response 
times. The time for brake torque to rise from 0 to the limiting value of 
15y400 foot-pounds (l8;l67.72 m-N) was 0.1 second for all three brakes and the 
value chosen for jimax runs. Also for all three brakes 

the acceleration value that caused a brake release signal to be generated was 
-50 radians per second per second^ and for brake application, zero radians per 
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second per second. For each brake; however; the time required for the torque 
to decrease from its maximum value to zero was different j that is, for brake A; 
t2 = 0.1; for brake B; t2 = 0.5; and for brake C; t2 = O.O5. The lowest 
value of brake decay time was used with brake C and is fairly representative of 
the fastest rate of brake release obtainable with current brakes. Although 
data on brake-response time are very scarce; it appears that modern brakes do 
not operate at frequencies much greater than 10 cycles per second. 

It can be seen in figure 8 that brake A reached a peak efficiency of 
92 percent when operated by a control having a frequency of about 5 cycles per 
second. The lowest efficiency for brake A of about 67 percent occurred at the 
lowest control frequency used for these calculations; that iS; 0.1 cycle per 
second. Brake B shows a much smaller variation in efficiency throu^out the 
control frequency range but its maximum value is significantly lower than that 
of brake A. The maximum efficiency of 49 percent for brake B occurs at the 
highest control frequency which was 100 cycles per second. Brake C exhibits 
the highest efficiency of all three brakes in the lower control frequency range j 
its efficiency rises from a "{l-'pevce'nt value at a control frequency of 0.1 cycle 
per second to about 90 percent at 0.5 cycle per second. When; however; the 
control frequency is increased beyond 0.5 cycle per second; the efficiency 
drops off very rapidly to a negligible value and indicates a complete loss in 
stopping ability. The reasons for the behavior of brakes A; B; and C; as shown 
in figure 8; can be explained by the computer time-history records obtained for 
some of the runs. 

Figure 9 shows the time histories obtained for brake A with an 
accelerometer-control frequency of 5 cycles per second; the condition of maxi- 
nrum efficiency for this brake. At the start of the braking cycle; it can be 
seen that the control acceleration values 3 are oscillating very rapidly and 
the frequency of this oscillation is about 90 cycles per second. Furthermore; 
the brake torque time-history trace indicates that the control must be gen- 
erating brake application and release signals since the brake torque never 
reaches its limited value. Therefore; the values of p must cover the range 
between 0 and -30 radians per second per second. Since the ratio of the fre- 
quency at which the accelerometer exercises control to its natural frequency 
is about 18 and thus would result in large attenuation to its response; it fol- 
lows that the wheel accelerations 0 which drive the accelerometer must be 
extremely large. This high-frequency control response enables the relatively 
low-frequency accelerometer to exercise precise; or very high response; control 
of the braking cycle; as will be shown shortly. It should be pointed out; how- 
ever ; that this response is dependent on the shape of the time variation in 
brake torque j that iS; if a more gradual or curved characteristic of torque 
plotted against time at the instant of application and release of brake torque 
had been used; the resulting values of 0 might not have been high enough to 
drive the accelerometer to the amplitudes required for brake control. Since 
the control is set to apply brake torque at zero angular acceleration and 
release it at -50 radians per second per second; brake torque is being applied 
over a larger part of each control cycle; therefore; the brake torque increases 
at the start of the run until the tire is operating at the maximum friction 
coefficient. From this point on; the tire begins to slip but since; as was 
just mentioned; the control is capable of very high response; this condition 
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prevents immediate lock-up of the wheel and, as indicated by the Sij time- 
history curve, controls the braking smoothly even though the tire is slipping, 

At the small slip ratios, control is well within the capabilities of the braking 
system since the slope of the p-slip curve for the tire is also small. As the 
slip ratio increases, however, the slope of the p-slip curve for the tire 
increases and tends toward a more unstable condition. (See fig. 5(a).) It 
can be seen that the accelerometer loses its ability to control the system 
smoothly when the tire slip ratio reaches a value of around 0.5- At this point, 
the entire character of the time histories are changed and the wheel and tire 
spin back up, and the cycle is repeated until a low rolling speed is reached. 

The decrease in efficiency for brake A, shown in figure 8, at the higher 
control frequencies is caused by the higher frequency accelerometer being cap- 
able of exercising smooth control over a larger portion of the back side of the 
p-slip curve. Therefore, the tire operates at higher slip ratios and at cor- 
respondingly lower coefficients of friction for a somewhat larger part of the 
run before instability occurs. This effect can be seen in figure 10 which shows 
the time histories obtained for brake A with a control accelerometer frequency 
of 50 cycles per second. 

The lower values of efficiency obtained for brake A at the lower control 
frequencies is basically caused by the low response of the accelerometer. It 
can be seen in figure 6(a), which is the computer run made by using brake A 
and a control frequency of O.5 cycle per second, that the time between brake- 
release signals, as indicated by the p trace and Tg trace, is about 0.2 sec- 
ond. This comparatively long time period allows the brake to lock the wheel 
and therefore causes the tire to operate at the full skid coefficient Psi^id 
throu^out a large proportion of the run. Since Pskid is a relatively low- 
friction value, the efficiency of the braking system is degraded. Although the 
motion p of the control is very complex and depends on the coupled response 
of the various masses and elastic properties of the system, a coarse physical 
understanding of the underlying reason for the very low response of the 
0‘5-cycle-per-second accelerometer as compared with the 5-cycle-per-second con- 
trol can be obtained from figure 11, Figure 11 shows the variation of p - 0 
as a f\inction of control natural frequency for the case of simple harmonic 
motion in which the maximum acceleration is 50 radians per second per second. 

It can be seen that for a natural frequency of O.5 cycle per second, the angu- 
lar displacement must be around 170° whereas for a natural frequency of 5 cycles 
per second, the excursion is less than 2°. 

The consistently low values of efficiency (between percent and 50 per- 
cent) exhibited by brake B throughout the range of control frequencies are due 
primarily to the relatively long brake release time t2 = 0.5 second. Fig- 
ure 12 shows the time histories of the run made with this brake by using a 
control frequency of 5 cycles per second. At the start of braking, the con- 
trol causes the system to reach Pmax I'apidly for much the same reason as 
described for brake A operating at this same control frequency. The rise in 
torque for this brake is actually faster than that for the equivalent brake A. 
The reason, of course, being that brake torque decayed at a much slower rate 
following the generation of a release signal. This slow rate of brake decay 


26 







Figure 11.- Variation of angular displacement^ of sensor control mass with respect 
to wheel, with natural frequency of sensor control for the condition of simple 
harmonic motion and maximum sensor angular acceleration of ^0 radians/second^. 

also prevented the accelerometer from exercising precise control at values of 
Srp > 0; therefore^ the wheel reached a full-skid condition in a relatively 
short period of time as compared with that for the eq^ui valent brake A. Fur- 
thermore^ the tire was not able to return to the nonslip condition after 
reaching full skid because of the slow brake -torque decay rate^ and thus oper- 
ated at or near jiski^ throu^out the remainder of the run. It should be 

pointed out that for this portion of the run^ the brake release signals were 
probably generated by the low wheel speed control. 

The very different variation of brake C with control frequency as com- 
pared with the other braking systems is basically due to the brake decay time 
t2 being significantly less than the application time tq. Figure 13 shows 
the run made with this brake at a control frequency of 0.2 cycle per second. 
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The efficiency, as indicated in figure 8, was about 79 percent, "^is braking 
system operates around the peak of the p,-slip curve for the tire (fig. 3(3’) 
and for most of the run, the excursions on both sides of this optimum point 
persist for only a relatively short period of time. This fact can be noted on 
both the &r and u time histories. The time history indicates that with 
the exception of the last several hundred feet of the run, the friction coef- 
ficient was held at a value which averaged close to hmax* This condition was 
made possible because, despite the low response of the control, the wheel was 
capable of driving the control to sufficiently high values of acceleration to 
generate brake release and application signals about 10 times each second 
during most of the run. Furthermore, as can be seen in figure 15, the response 
of the control signal appeared as a distorted wave and its shape was such that 
brake torque was being applied during a greater portion of the time of each 
braking cycle. Despite the longer time of brake application, the tire was not 
able to reach a large value of slip ratio because of the rapid rate at which 
brake torque decreased once a release signal was generated. The highest effi- 
ciency obtained for brake C occurred at a control frequency of 0.5 cycle per 
second and was a result of this control generating more brake release signals 
per unit time than for the brake operating at 0.2 cycle per second; and thus 
the excursions of the friction coefficient on each side of the maximum were 
shortened. This effect can be seen in figure 1^4-. 

The very abrupt drop in braking efficiency that occurred when the control 
frequency was increased beyond 0.5 cycle per second again resulted from the 
fact that the amplitude response of this control to the wheel angular accelera- 
tions 0 was sufficient to generate brake release signals at a very high fre- 
quency. Furthermore, these high-frequency control accelerations appeared as a 
fairly pure wave form, and since the brake torque decayed at double the rate 
at which it was applied, it was not possible for this system to develop any 
brake torque. This effect can be seen in figure 15 which shows only a short 
section of the run made with brake C at a control frequency of 10 cycles per 
second. 

The foregoing discussion of figure 8 indicates that the efficiency of a 
braking system depends on a rather complex interaction between the values of 
tire frequency, brake response, control frequency, and the slope of the p-slip 
curve for the tire. Since the combinations investigated indicated that the 
efficiency of braking systems could vary from more than 90 percent down to 0, 
this result would Indicate the need for careful analysis and testing in the 
selection and integration of the individual hardware units which comprise the 
braking system. These results also indicate that under certain dynamic condi- 
tions, an accelerometer with a relatively low response can exercise precise 
control of brake torque over a large part of the back side of the ^x-slip curve 
for the tire. 


Low Tire-Runway Friction Coefficients 

Causes of low friction coefficients .- As was indicated previously, current 
braking systems operate satisfactorily on dry or even damp runway surfaces. 

For these cases, substantial retarding forces can be developed by the braking 
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second; t 2 = 0.05 second; — ^ = 0.5 cycle/second; fimax 






Figure I5,- Analog computer time-history solution using brake C. = Oj 

p’2 = -30 radians/second^ j t]L = second; t2 = O.O5 second; 

= 10 cycles/ second; = O.U. 


system throughout the landing run. However ^ braking systems must also operate 
when the tire -runway coefficient is very low; the most common cause for low 
friction coefficients is the presence of water or slush on the runways. For 
such conditions^ a phenomenon commonly called tire hydroplaning or aquaplaning 
can occur. During complete hydroplaning, the airplane tire is completely sup- 
ported by the fluid and the wheels may come to a complete stop even though the 
airplane has an appreciable forward velocity and no braking torque is applied. 
Much information regarding this phenomenon can be found in the literature. 

(See, for example, refs. 6, 8, and 9-) 

In brief, complete hydroplaning occurs at a velocity which depends pri- 
marily on the tire pressure, tread pattern, and fluid depth. For a modern 
large commercial jet transport having a tire inflation pressure of I50 Ibf/in 
(103.4 N/cm^), the velocity at which complete hydroplaning begins is approx- 
imately 110 knots, and for the tire used in this analysis complete hydroplaning 
occurs at about I50 knots. As the forward speed drops below the critical speed 
for complete hydroplaning, the wheel once again contacts the runway surface, 
and the load supporting the tire is then divided between the fluid and the run- 
way. The distribution of load between the fluid and runway for this partial 
hydroplaning condition depends primarily on the magnitude of the forward speed. 
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Figure 15.- Concluded. 


The inanner in which the load is distributed between the fluid and runway deter- 
mines the amount of retarding force that the braking system can develop; that 
is^ at high forward speeds where a small portion of the load is carried by the 
runway^ the braking forces that can be developed are smaller than those that 
can be developed at the lower speeds where more of the airplane weight is sup- 
ported by the runway. Under forward velocity conditions resulting in complete 
hydroplaning^ it is clear that braking systems are completely ineffective and^ 
as a result^ during this time the aerodynamic drag developed by the airplane 
spoilers and wing flaps and reverse engine thrust must be depended on to reduce 
the speed. Once the speed has been reduced to a level that allows the tires to 
make contact with the runway^ however^ the braking system may then become 
effective in two ways in bringing the airplane to a safe stop: first, by devel- 

a sufficient retarding force throughout the entire distance to stop the 
airplane in the available runway length with an adeq[uate safety margin; second, 
by allowing the tires to operate at a minimum slip ratio so that cornering 
forces can be developed by the rolling tires to maintain directional stability 
which is lost when the wheels lock up or when the tires operate at high slip 
ratios. Hydroplaning, as pointed out, depends on velocity; however, other 
causes of low tire-surface friction coefficients, such as those for icy sur- 
faces, are, in general, quite independent of velocity and exhibit a low coef- 
ficient throughout the entire velocity range. 
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Braking; sysiein efficiancy.- In order to study the effect of variations in 
on the efficiency and performance of braking systems, a series of com- 
puter runs were made over a range of varying from 0.1 to 0.8. 

Since the previous results Indicated that brake A operating with a control fre- 
quency of 5 cycles per second was the most efficient system for M-max ~ 0.4 
(see fig. 8), this system was used for the first series of runs made over the 
coefficient-of -friction range. Results showing the efficiency of this braking 
system are presented in figure l6 by the curve drawn throu^ the circle symbols. 
For the runs made at M-max values greater than 0.5^ the limiting brake torque 
of 15,400 foot-pounds (18,167.72 m-N) is less than the product there- 

fore, during these runs the system operated at about the maximum brake torque 
throughout the entire run with no tire slipping, and the stopping distances 
were about the same. The apparent Increase in efficiency associated with the 
decrease in pj^ax caused by an increase in the value of 

minimum runout distance used in computing the efficiency parameter. As Pmax 
decreases below 0.4, it can be seen that the braking- system efficiency also 
decreases in this range and reaches a value of less than 50 percent when 
p =0.1. This effect is Just opposite the effect required to Insure safe 

aircraft operations. To illustrate this point, the efficiency required to stop 
an airplane over the given coefficient of friction range in 'J ^000 feet 
(2,155.60 m) is shown in figure I6 as a dashed line. For Pmax ~ 0.8, an effi- 
ciency of less than 20 percent is required to stop in 7,000 feet (2155.60 m) 
whereas for Pj^^x = efficiency of about 90 percent is required. The 

stopping distance computed for the run (represented by the circle symbol) at 



_l -j -L 1 J X X- 

.1 .2 ,3 .4 .5 .6 .7 


Maximum coefficient of friction, 


Figure l6.- Effect of varying iDrake-control signal on the 
variation of braking- system efficiency with maximum fric- 
tion coefficient, tx = t 2 - 0.1 second; 

^ 5 cycles/second. 
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lajjiax = 0.1 was around 15,000 feet (4,572 m) . This result indicates that a 
braking system which operated satisfactorily in the presence of relatively high 
or medium coefficients of friction could operate unsatisfactorily at low coef- 
ficients of friction. 

In seeking a corrective measure which would involve minimum modifications 
to the braking system, a series of computer runs were made over the range of 
friction coefficients with the same brake but changing the value of the 

acceleration at which the control generated a brake-release signal^.^ For these 
runs P 2 made equal to -5 radians per second per second and was the 

same as that for the previous runs, that is, zero. The results of these com- 
puter runs are indicated by the square symbols in figure I 6 . This change 
resulted in essentially a runout of somewhat less than 7^CK30 feet ( 2 , 153.60 m) 
for each run over the entire range of computer time 

histories for both runs appearing in figure I 6 that were made at ^inax = 0.2 
are shown in figures I 7 and I 8 . 

Figure I 7 shows the run made with P 2 = "50 radians per second per second. 
For the first several thousand feet of runout, braking torque is being cycled 
very rapidly by the control in the manner di earlier. The control is 

again capable of operating smoothly on the back side of the p.-slip curve for the 
tire and continues to exercise control until the tire is, for all intents and 
purposes, in a full skid. For this low condition, however, the wheel 

is not able to spin up to full rolling speed and the control continues to 
operate at this point for the remainder of the run which is very close to 
Srp = 1. This condition, of course, results in much tire skidding. Again, the 
brake-release signal is probably generated by the low wheel speed control. 

In contrast, figure I 8 shows that the run made with 
P 2 = -5 radians per second per second operates throughout the entire stopping 
distance with no tire slipping. The fact that for both runs the runout dis- 
tance or efficiency was practically the same indicates that both operated over 
the same average value of friction coefficient. The difference being that for 
the run in which P 2 = -30 radians per second per second, the braking system 
operated on the back side of the fi-slip curve for the tire at |i ~ l^skld^ 
whereas although the other system operated at about the same average friction 
value, its operation was on the front side of the (i-slip curve for the tire 
where no tire skidding occurs. 

This optimum type of operation can be explained by noting that the rolling 
acceleration of the wheel at -5 radians per second per second is equivalent to 
a linear forward acceleration of -0.l8g. Since djnax 0.2, it is evident 

that brake torque was released just before the tire reached the limiting fric- 
tion value for the runway and although the high-frequency oscillations of 3 
are barely perceptible on the record shown in figure I 8 because the value is 
so small, it is evident that the brake torque was controlled at a value that 
kept the tire rolling without slipping at an apparent coefficient of friction 
just under M-niax* other words, having P 2 = "5 radians per second per second 
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resulted in a linear deceleration x slightly under that specified hy equa- 
tion (l) for a |i value of 0,2. 
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Figure l8.- Analog computer time-history solution. 3]^ = 0} P 2 = "5 radians/second^; 

C% 

%x - t 2 = 0.1 second; = 5 cycles/second; = 0.2. 


Figure I9. shows the run made at (iniax ~ with the braking system set 

so that 32 “5 radians per second per second. Although the tire is slipping 

throughout most of the run because M-j^ax less than 0,l8, the increase in 
slip is very gradual^ and as indicated in figure l6^ the efficiency for this 
system is over 90 percent. 
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Braking- system performance . - The effectiveness of the two hraklng systems 
In minimizing tire skidding over the range of friction coefficients is indicated 
in figure 20. For the braking system in which 32 set at -50 radians per 


second per second, the percentage of stopping energy contributed by tire 
skidding rises sharply as the available coefficient of friction decreases , and 

where l-iiiiax ~ 0*1^ almost 



Maximum coefficient of friction, 


Figure 20.- Effect of varying brake-control signal on 
the variation of work done by tire in skidding with 
maximum friction coefficient, tj = t2 = 0.1 second; 


80 percent of the stopping 
energy is contributed by 
tire skidding. On tbe 
other hand^ for the system 
where 02 = "5 radians per 
second per second; tire 
skidding is evident only 
for the minimum coefficient 
of friction of 0-1 and for 
this case contributes only 
slightly more than 10 per- 
cent of the total stopping 
energy. 

When the performance 
of a braking system is 
evaluated with respect to 
its ability to limit tire 
skidding, it is iniportant 
to associate the perform- 
ance with the magnitude of 
the available coefficient 


— ^ = 5 cycles/second. 
2it 


of friction. Skidding at 
high coefficients of fric- 
tion results in excessive 

tire wear; whereas diiring skidding at low coefficients of friction, the tire 
wear is greatly reduced. The importance of maintaining skidding at a minimum 
during stops at a low coefficient of friction condition is associated with 
the loss of direction control. During skidding the available cornering force 
is decreased and becomes zero for a fully locked tire. Needless to say, this 
condition is a dangerous one, especially if complete locking occurs early in 
the landing when the forward speeds are greatest. The hazards associated with 
this condition are discussed in reference 10. 


Modulation of brake-release signal .- The results shown in figures l6 and 
20 would suggest that a hraklng system which modulated the value of 32 ^ 

fiinctlon of airplane forward speed could be very effective in achieving safe 
stops on flooded or slush-covered runways. Such a system could be adjusted so 
that the deceleration developed by the flaps and spoilers in addition to engine 
reverse thrust was greater than the initial value of 32 that would cause it 
to release brake torque. As a consequence, no hraklng would he applied during 
the early stages of the landing since the wheel deceleration would he greater 
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than ^ 2 ' Therefore^ the full cornering capability of the freely rolling wheels 
would be available to maintain directional control of the aircraft during this 
critical stage of the landing runout. Reference 11 indicates that althou^ the 
aerodynamic drag becomes less as the forward speed decreases^ the average 
deceleration during the early stages of the landing runout attributable to both 
aerodynamic drag and engine reverse thrust is around 0.2g. As the velocity 
decreases and the amount of the airplane weight supported by the runway 
increases, the stopping effort contributed by the brakes could then be augmented 
by modulating the value of P2 to allow the release of brake torque at a larger 
wheel deceleration. 

For the conditions of low friction coefficient which are not largely 
dependent on velocity, such as ice, the foregoing modulation system could not 
be used since ^^uld be practically the same for the entire braking stop. 

It would then be necessary to maintain the original small value of P2 ^*or the 
entire run. The braking process, however, could still be accomplished with 
high efficiency, as indicated by the curve associated with the square symbols 
in figure I6. Dispensing with the modulation of ^2 wuld, of course, affect 
the stopping distance. 

In this connection, it should be mentioned that low coefficients, which 
are also largely independent of forward velocity, can be experienced on very 
smooth damp surfaces. In this case, the low coefficient results from a lubri- 
cation process. It was shown in reference 12, however, that by merely making 
such surfaces rough, very large increases in the coefficient -of -friction values 
can be realized. It would, therefore, appear that the braking- system problem 
could be substantially alleviated by insuring that runway surfaces have ade- 
quate roughness. 

Efficiency . - It should be pointed out that the method which has been used 
in defining efficiency might be difficult to use in practice. The reason for 
this difficulty lies in the fact that the numerator of the efficiency term, 
minimum stopping distance, must be associated with a definite coefficient of 
friction. Therefore, efficiency, as defined herein, would only take on real 
meaning when defined as a curve over a specific range of coefficients of fric- 
tion. The dashed curve in figure I6 is an example of one type of efficiency 
curve and would be the efficiency curve for a constant stopping distance of 
7,000 feet (2,135.60 m). 

Variation of p with Sji for Sfp > 0.- It was pointed out in the dis- 
cussion concerning operation of brake A that the slope of the p-slip curve for 
the tire in the region where Srp > 0 affected the stability of the braking 
system and therefore had a significant effect on the efficiency and operation 
of the system. Also, the actual variation of p with Sip for the region in 
which Srp > 0 could not be described by a unique curve because it depended on 
a number of variables, some of which varied in a random manner. The slope of 
the curve used to obtain all the results presented was based on an empirical 
value determined from experimental data obtained during dry-surface braking 
tests, and the value of about one-third of the dj^^x 
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should he pointed out, however, that the available experimental data appear to 
indicate that when braking stops are made on surfaces which exhibit low coeffi- 
cients of friction, the values of M'skid ^^max differ in value as 

much as the dry-surface values. This difference would, of course, affect the 
slope of the ^-slip curve for the tire. This result would suggest the need for 
further studies to establish critical variations of la with St compatible 
with the conditions encountered in practice for use in the design of braking 
systems . 


CONCLUDING REMARKS 


A preliminary analog computer study has been made of the operation of an 
automatically controlled braking system in which the elastic character of the 
tire and the response time of the brake were taken into account in developing 
the system equations of motion. These equations were derived in general terms 
for a simplified single wheel and tire system and then solved on the analog 
computer for a particular case of a wheel equipped with a 32 X 8.8 type VIl 
tire. A relation was developed to define the variation of the coefficient of 
friction with slip ratio and this relation was in agreement with some experi- 
mental results. The brake simulation used was torque limited, and the varia- 
tion of brake torque with time during application and release was assumed to 
be linear. The solutions indicated that tire elasticity has a very important 
effect on the operation of the system. For the case of an elastic tire, the 
wheel motion could differ considerably from the tire motion at various instants 
of time during the braking process. As a result, during cyclic braking, the 
variation of the coefficient of friction with wheel slip ratio could have both 
positive and negative slopes in the region between free rolling and the attain- 
ment of maximum friction coefficient. This result is in contrast to current 
notions since experimental data are normally presented with this slope drawn 
positive. The results also indicated that the efficiency or stopping capa- 
bility of an automatic braking system is affected by tire elasticity. 

Studies were also made of several brakes having a different variation of 
torque decay with time. It was found that the braking- system efficiency 
depended on a complex interaction between the values of tire frequency, brake- 
response time, and skid-control sensor frequency. The results showed that 
depending on the combinations of the foregoing factors, the efficiency of the 
braking systems analyzed varied from better than 90 percent down to 0 percent. 

It was also found that under certain dynamic conditions, a skid-control sensor 
with a relatively low response was capable of exercising control of brake torque 
at a frequency about l8 times its own natural frequency. 

Studies were also made of the effect of low tire-ground friction coeffi- 
cients on the operation of automatic braking systems. The results showed that 
a braking system which operated satisfactorily in the presence of relatively 
high or medium coefficients of friction could operate unsatisfactorily at low 
coefficients of friction. For low friction coefficients, large decreases in 
both stopping distance and tire skidding could be realized by proper selection 
of the wheel deceleration value used to generate a brake-release signal. The 
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results suggested that modulation of this control signal could he very effec- 
tive in achieving safe stops on water or slush -covered runways. 

The results of this preliminary study cannot^ however, he considered as a 
precise description of control-system operation since many factors affecting 
the operation of a braking system are not presently known. Although the trends 
indicated herein can he considered valid, further research is needed in areas 
related to the general mechanical properties of rolling tires, the character of 
the loads developed hy such tires under the influence of brake torque, and the 
time characteristics of brake-torque variation. A better understanding of these 
phenomena would result in Improved accuracy of the analog simulation and the 
results would then be more descriptive of those obtained during actual braking 
stops. The results also suggested the need for further studies to establish 
critical variations of coefficient of friction with tire slip ratio that would 
be compatible with conditions encountered in practice. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 
Va., August l6. 


1965. 
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APPENDIX A 


DERIVATION OF EQUATIONS OF MOTION 


The equations of motion for the mathematical model shown in figure 1 will 
be obtained using La Grange's equations which are (ref. I 3 ): 


dt \dq^ y Sqj^ 

The kinetic energy of the system is 


(1 = x,9jCp,3) (Al) 


T 


2 g 



•P "I *2 

2 ^T^ ^ 



The potential energy of the elastic elements of the 


system are 


(A2) 



The generalized forces are 


(A3) 


(aA) 


Qx = -hW (A5) 

^0 ~ ^ (A6) 

= ripjiW — r^jCq)(cp - 0) (A7) 

Qp = -nCp(p - 0) (A8) 


The sign of Tg is chosen so that the brake torque always opposes the wheel 
motion 9 . 


k& 



APPENDIX A 


The specific equations of motion when 1 = x, 0, cp, and p are as 
follows : 

When i = x: Substituting equations (A2) and (A5) into equation (Al) gives 


X = -pg 


(A9) 


When 1=0: Substituting equations (A2), (Aj), (a 4), and (a 6) into equa- 

tion (A1) and noting that 


= ^w^T - 9) 

and 

yp = n(p - 0) 

gives 


1^0 - l^r2(cp - 0) - kpn2(p - 0) = Tg + t ^ q {^ - 5) 


Since the torque produced by the deflection of the control spring is negligible 
compared with the torque produced by the elastic tire^ the term -Icpn^Cp - 0) 

is dropped and this equation can be written as 


1^0 - krpr§(cp - 0) 


I’d ’’w^0(t - 0) 


Also 


(AlO) 


Cc,0 


2 ^ W*^,0 


(All) 


where 



Dividing both sides of equation (AlO) by equation (All) and then multiplying 
both sides by 2co^^0/r^ gives 
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APPENDIX A 


T 

0 - 2o)n,0le((p - 0) - o:g,Q(cp - 9) = (A12) 

% 

where 

!e = (A13) 

Cc,9 

When i = cp: Substituting equations (A2), (A5)^ and (A7) into equa- 

tion (A1) gives 


Irpcp + torr^(cp - 0) = r^W - r^cp(cp - 0) 


(A14) 


Also 


where 


^c,0 




(A15) 


«4i,cp = 



Dividing both sides of equation (A14) by equation (AI5) and then multiplying 
both sides by gives 


9 + 2c%,cplcp('P - S) + - 9) = 


It 


w 


(A16) 


where 


Icp = 



By noting that C^p = Cq , it can be shown that 


^ ~ ^0 



(AI 7 ) 


(AI8) 
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When i = p: Substituting equations (A2), (a 4), and (A8) into equa- 

tion (A1) gives 


Ipp + k|3n2(p - e) = -nCp(p - 0) 


(A19) 


Also 


^C,p 




g. P-ui 

n 


(A20) 


where 




kc,n^ 


Dividing both sides of equation (AI 9 ) by equation (A20) and then multiplying 
both sides by 2o)jj^p/n gives 


where 


p = -2a^^p^(p 


0) - o:g^p(p - 0) 



(A21) 


(A22) 


Equations (A 9 ) j, (A12), (AI 6 ), and (A2l) are the equations of motion of 
the system and are the same as equations (l), (2), (3)^ aiid. (4). 
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APPENDIX B 


DERIVATION OF BRAKE -TORQUE EQUATION FOR THE CASE OF A FULLY K)CKED WHEEL 


The brake-torque equation for the case of a fully locked wheel is derived 
in this appendix. 


During complete wheel locking 


0=0 = 0 


Also, from appendix A 


and 


2 

“n,0 


% 



It 


(Bl) 


(B2) 


(B5) 


Substituting equations (Bl), (B2), and (BJ) into equations (A12) and (A16) and 
then adding the two resulting equations gives 

+ 2cp^o^^q)lTlcp - “n,0lW?0) “ U^^T = 

where Tg is the symbol for the value of brake torque during the time the 
wheel is fully locked. 

By using the expressions of appendix A for <Pn,cp^ ^,0^ ^0^ 

it can be shown that the coefficient of 2cp (the expression in parenthesis) 
in equation (Bif) is equal to zero. Therefore, the expression for evaluating 
brake torque while the wheel is locked is 

Tg^g = It<P - M-Wrtp (B5) 

Equation (B5) is identical to equation (8). 
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APPENDIX C 


DERIVATION OF THE EQUATION FOR DETERMINING THE WORK DONE 
BY THE TIRE IN SKIDDING 

THg work d.on6 In skidding is 6Q^u.3.1 to ill© produ.ci of ill© skidding disisnc© 
and tH© ground drag force. Expressed in incremental form: 


^^skid = ^skidl^'^ 

Hut 

^skid = 

Substituting equation (C2) into equation (Cl) gives 

^^skid. ” xSij^WAt (^3) 

Dividing botb sides of* equation (Cj) by At and letting At approach zero as 
a limits that is 


limit = iSTnW 

At-^0 


Therefore 




skid 


pt=t 

= w / 

^ t=0 


xSip dt 


Equation (c4) is identical to equation ( 9 ). 


(C4) 
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